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LONG-TERM  GOALS 

The  long  term  goal  is  to  develop  a  robust  multi-phase,  multi-elass  numerieal  modeling  framework  for 
both  cohesive  and  non-cohesive  sediment  transport  in  the  fluvial  and  coastal  environments. 

OBJECTIVES 

The  objectives  of  the  present  study  focus  on  extending  a  fluid  mud  model  for  boundary-layer-driven 
and  gravity-driven  transport  (Hsu  et  al.  2007)  with  several  new  capabilities.  Specifically,  the  objectives 
are  to: 


•  Conduct  comprehensive  model-data  comparisons  with  laboratory/field  observed  fluid  mud 
processes,  specifically  in  conjunction  with  new  MURI  initiative  for  understanding  wave-mud 
interaction. 

•  In  conjunction  to  NOPP-Community  Sediment  Transport  Model  initiative,  the  fluid  model 
is  utilized  to  provide  parameterizations  for  wave-boundary-layer-scale  transport  processes. 

•  Extend  the  existing  fluid  mud  model  with  a  bed  module  to  enable  direct  modeling  on 
consolidation/fiuidization.  Extend  the  model  to  multi-dimensional  for  various  coastal  and 
estuarine  applications. 

APPROACH 

It  is  important  to  understand  the  fate  of  terrestrial  sediment  into  the  coastal  ocean,  because  it 
determines  for  example  the  seabed  properties  and  the  turbidity  of  the  water  column.  Several  recent 
initiatives,  such  as  NOPP  Community  Sediment  Transport  Model,  MURI-  Understading  Wave-Mud 
Interaction  and  Tidal  Mud  Flats  DPI  have  put  forward  new  outstanding  science  and  technical 
questions  for  coastal  sediment.  The  success  of  these  new  studies  inevitably  depends  on  our  level  of 
understandings  on  the  crucial  small-scale  mechanisms  that  can  be  either  resolved  or  parameterized  via 
detailed  measurement  and  modeling.  These  new  initiatives  also  have  a  similar  goal  to  establish 
improved  understanding  on  coastal  hydrodynamics  and  sediment  transport  processes  in  heterogeneous 
environment  with  more  emphasis  on  cohesive  sediment  transport. 
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The  dynamics  of  fluid-mud  transport  involve  a  variety  of  physical  mechanisms,  including  for  example, 
the  boundary  layer  and  gravity-driven  transport,  turbulence  modulation,  flocculation,  non-Newtonian 
rheological  behavior  and  consolidation  (Dyer  1989;  Mehta  1989).  Hence,  a  general  modeling 
framework  for  fluid  mud  need  to  be  based  on  multiphase  flow  theory.  In  this  study,  a  fluid  mud  model 
based  on  Fast  Equilibrium  Eulerian  Approximation  (Eerry  and  Balachandar  2001)  to  the  multiphase 
equations  appropriate  for  fine  sediment  has  been  developed  and  extended  to  model  various  cohesive 
sediment  transport  processes.  An  earlier  version  of  this  fluid  mud  model  (Hsu  et  al.  2007)  has  been 
shown  to  be  capable  of  modeling  tidal-driven  fluid  mud  (Kineke  et  al.  1996)  and  wave-supported 
gravity-driven  mudflows  (Traykovski  et  al.  2007). 

WORK  COMPLETED 

Eurther  development  of  the  fluid  mud  model  and  understanding  of  various  cohesive  sediment  transport 
processes  are  carried  out  in  several  directions: 

Extend  to  2D  and  incorporate  Bingham  rheology:  The  previous  IDV  fluid  mud  modeling 
framework  (Hsu  et  al.  2007),  is  extended  to  2D  to  simulate  wave-mud  interaction.  A  2D  numerical 
wave  model  (COBRAS,  Ein  &  Eiu  1998)  based  on  Reynolds-averaged  Navier-Stokes  (RANS) 
equations  with  Volume  of  Eluid  (VOE)  method  for  free-surface  tracking  is  modified  to  solve  the 
proposed  fluid  mud  equations.  Notice  that  the  proposed  governing  equations  for  fluid  mud  and  the 
related  closures  reduce  to  single-phase  RANS  equation  when  mud  concentration  becomes  zero.  Hence, 
the  numerical  model  calculates  wave  propagation,  sediment-laden  wave  boundary  layer,  fluid  mud 
transport  and  their  interactions  continuously  and  consistently  with  one  single  set  of  balance  equations 
and  closures  (i.e.,  no  matching  is  required  at  mud-water  interface).  Currently,  we  calculate  the 
rheological  stress  based  on  a  mathematical  form  that  is  similar  to  Bingham  rheology  (Ee  Hir  et  al. 
2001;  Mei  and  Eiu  1987;  Erigaard  and  Nouar  2005).  Essentially,  a  sediment  viscosity  is  utilized  which 
is  depending  on  local  strain  rate  and  fluid  mud  concentration.  Model  results  for  cnoidal  wave 
propagating  over  a  190  meter  muddy  seabed  are  shown  in  figure  1. 

Parameterization  of  wave-supported  gravity-driven  mudflows:  Typical  coastal  models,  such  as  the 
ongoing  NOPP-Community  Sediment  Transport  Model  (NOPP-CSTM)  are  not  designed  to  resolve  the 
thin  wave  boundary  layer  near  the  seabed.  Hence,  processes  that  occur  within  the  wave  boundary  layer 
need  to  be  parameterized  as  sub-module  to  provide  quantities  such  as  bottom  drag  coefficient  and 
sediment  transport  rate.  The  fluid  mud  model  is  utilized  to  study  fluid  mud  transport  in  the  wave- 
current  boundary  layer.  Specifically,  model  results  are  used  to  provide  parameterizations  for  wave- 
supported  gravity-driven  mudflows  (Traykovski  et  al.  2007).  The  nature  and  the  existence  of  wave- 
supported  gravity-driven  mudflows  are  diagnosed  by  varying  the  floe  properties,  bed  erodibility  and 
rheological  stresses.  By  selecting  representative  downslope  fluid  mud  transport  events,  numerical 
experiments  are  conducted  to  study  the  effect  of  wave  and  long-shelf  current  intensities  on  the  fluid 
mud  transport,  and  its  parameterization  (Wright  et  al.  2001).  Model  results  suggest  that  the  drag 
coefficient  decreases  with  increasing  wave  intensity,  and  seems  to  follow  a  power  law.  The  bulk 
Richardson  numbers  is  less  sensitive  to  wave  intensity  but  has  a  magnitude  about  0.08,  which  is  factor 
3  smaller  than  the  critical  Richardson  number.  The  dependence  of  wave-supported  gravity-driven 
mudflows  on  long-shelf  current  intensity  is  also  less  sensitive  compared  to  that  on  wave  intensity. 
However,  when  the  long-shelf  current  is  as  large  as  about  1.0  m/s  at  1  m  above  the  bed,  auto¬ 
suspension  may  occur. 
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Flocculation:  A  balance  equation  for  floe  diameter  in  response  to  floe  aggregation  and  break-up 
proeess  in  turbulent  flow  (Winterwerp  1998)  is  ineorporated  in  the  IDV  fluid  mud  model.  The 
numerieal  model  is  able  to  ealeulate  the  spatial  and  temporal  variation  of  floe  size  in  turbulent 
boundary  layer  (with  turbulenee  quantities  provided  by  the  k-s  elosure)  using  a  fixed  fraetal  dimension. 
Notiee  that  fraetal  dimension  is  required  to  speeify  the  exeess  density  of  the  mud  floe  (and  henee  the 
settling  veloeity).  Based  on  field/laboratory  observations,  the  fractal  dimension  is  not  a  eonstant  and 
shall  ehange  dynamieally  with  the  earner  turbulent  flow  (e.g.,  Dyer  and  Manning  1999).  Reeently, 
Khelifa  and  Hill  (2006)  proposed  a  model  for  settling  veloeity  that  utilizes  a  variable  fraetal  dimension 
depending  on  the  floe  size  itself  (e.g.,  fraetal  dimension  beeomes  3.0  when  floe  size  reduees  to  that  of 
primary  partiele).  They  demonstrate  that  variable  fraetal  dimension  is  more  appropriate  to  relate 
measured  floe  size  with  settling  veloeity  in  the  field  and  laboratory.  In  this  study,  we  further  derive  the 
balance  equations  for  floe  dynamics  based  on  the  variable  fraetal  dimension  and  eome  up  a  new  model 
for  floe  dynamies  that  is  able  to  be  eoupled  with  numerieal  model  for  fluid  mud  transport.  The  new 
formulation  is  shown  to  prediet  the  equilibrium  floe  size  under  different  turbulent  shear  rates  measured 
in  the  laboratory  (Son  and  Hsu,  submitted). 

Development  of  a  unified  transport/consolidation/fluidization  model:  Traditionally  the 
geoteehnieal  engineering  eommunity  has  been  interested  in  the  eonsolidation  proeesses  using  for 
example  Gibson’s  equation  (Gibson  et  al.  1967).  However,  in  order  to  predict  morphological  evolution 
in  a  dynamic  muddy  environment,  both  the  consolidation  and  fluidization/liquefaction  under  turbulent 
shear  flow  have  to  be  modeled.  Henee,  there  is  a  need  to  eome  up  with  a  more  general  modeling 
framework  for  mud  bed  dynamies.  In  terms  of  the  governing  equation,  existing  simplified  multiphase 
flow  theory  ean  be  direetly  adopted  (see  Toorman  1999  for  a  similar  approaeh).  The  diffieulties  are  the 
appropriate  closures  for  rheology,  effeetive  stress  and  permeability.  A  series  of  work  by  Kranenburg 
(1994)  and  Merekelbaeh  and  Kranenburg  (2004)  suggested  that  if  we  assume  floe  is  self-similar  on  all 
seales,  speeifie  relations  exist  between  the  fraetal  dimension  and  various  properties  of  the  eohesive 
sediment,  sueh  as  the  rheology,  effeetive  stress  and  permeability.  Applying  sueh  theory  to  existing 
eonsolidation  test  data  Kranenburg  (1994)  found  that  fraetal  dimension  is  around  2.75,  whieh  is  mueh 
larger  than  that  of  typieal  floe  aggregates  in  mobile  suspension.  This  suggests  that  a  eomplete 
eonsolidation/fiuidization  formulation  shall  incorporate  floe  dynamies  with  variable  fraetal  dimension. 
We  are  eurrently  working  on  the  theoretieal  formulation  of  this  process.  In  the  future,  such  formulation 
can  be  direetly  adopted  in  the  fluid  mud  model. 

RESULTS 

The  fluid  mud  modeling  framework  is  implemented  into  a  numerieal  wave  model  COBRAS  to  direetly 
simulate  wave-mud  interaetion.  Here,  we  deseribe  a  simple  example  to  demonstrate  the  model 
eapability.  A  Cnoidal  wave  train  of  wave  height  0.72m,  wave  period  6  see  is  sent  into  the  numerieal 
wave  tank  of  water  depth  2.5  m  and  total  length  200  m  (see  Fig  1).  Between  x=10  m  and  190  m,  fluid 
mud  of  floe  size  d=22  mp  and  speeifie  gravity  5=1.34  is  allowed  to  be  suspended  by  waves.  A 
snapshot  of  the  model  results  at  t=12  see  is  shown  in  Fig  la  for  the  entire  numerieal  wave  tank  and  Fig 
lb  for  a  elose-up  view  near  the  bed.  A  thick  layer  of  fluid  mud  is  suspended  with  a  thiekness  of  about 
50  em.  Fine  grid  of  zlz=2mm  in  the  vertieal  direetion  is  implemented  and  henee  the  wave  boundary 
layer  strueture  and  lutoeline  are  resolved.  The  grid  size  in  the  horizontal  direetion  is  zh:=20om,  whieh 
is  appropriate  to  resolve  the  Reynolds-averaged  wave  field.  The  dashed  eurve  near  the  free-surfaee  in 
Fig  la  represents  the  free-surfaee  elevation  from  another  independent  eomputation  without  the  fluid 
mud  (i.e.,  elear  fluid  case)  at  the  same  instant.  Apparently,  the  wave  height  starts  to  deeay  as  the  wave 
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train  propagates  further  downstream  over  the  muddy  bottom.  Fig  le  presents  model  result  of  time 
history  of  free-surfaee  elevation  at  wave  guage  x=190  m  (solid  eurve).  This  result  ean  be  eompared 
with  guage  data  obtained  at  x=10m  (right  before  the  waves  propagate  into  the  mud  regime,  light-dotted 
curve)  or  obtained  at  the  same  location  (x=190  m)  but  due  to  another  independent  computation  without 
the  presence  of  mud.  We  observed  that  the  wave  height  is  damped  by  about  18%  due  to  propagation 
over  muddy  bottom  of  180  m  in  length. 

More  insights  can  be  revealed  when  examining  the  near-bed  wave  boundary  layer  structure  (see  Fig 
lb).  Without  the  fluid  mud,  the  thickness  of  the  wave  boundary  layer  is  only  about  2cm  (e.g.,  see 
dashed  curve  atx=23  m),  which  is  about  10  time  smaller  than  that  with  the  fluid  mud  (about  20cm,  see 
solid  curve  at  x=23m).  Due  to  the  presence  of  fluid  mud,  large  rheological  shear  stress  provides  an 
additional  mechanism  to  mix  the  flow  momentum  across  the  wave  boundary  layer  (in  additional  to 
turbulent  shear  stress)  and  hence  the  resulting  wave  boundary  layer  is  much  thicker.  This  enhanced 
wave  boundary  layer  further  allows  suspension  of  a  thicker  layer  of  fluid  mud  (note  that  the  thickness 
of  fluid  mud  layer  often  scales  with  wave  boundary  layer  thickness).  The  damping  of  the  wave 
amplitude  is  apparently  caused  by  (1)  direct  energy  dissipation  through  rheological  stress  term  in  the 
momentum  equation;  (2)  indirect  mechanism  due  to  greatly  enhanced  wave  boundary  layer  thickness 
(and  hence  turbulence  in  the  wave  boundary  layer).  In  summery,  the  damping  of  wave  due  to  fluid 
mud  is  rather  interactive.  The  thicknesses  of  the  wave  boundary  layer  and  fluid  mud  layer  are  inter¬ 
dependent  and  it  is  difficult  to  prescribe  it  as  a  priori. 

More  comprehensive  simulations  using  various  wave  conditions,  mud  properties,  rheological  closures 
shall  be  conducted  in  year  2.  Currently,  the  computation  is  limited  by  CPU-time.  In  year  2,  a  postdoc 
researcher.  Dr.  Torres-Freyermuth,  who  use  to  develop  a  more  efficient  version  of  COBRAS  (Dr. 
Losada’s  group  in  Universidad  de  Cantabria,  Spain)  based  on  Fortran90  with  C++  library  will  be 
joining  UF  and  work  on  the  2D  fluid  mud  model  for  wave-mud  interaction.  In  year  2,  we  will  also  start 
working  model-data  comparison  on  measured  wave  boundary  layer  and  fluid  mud  data  by  co-PI 
Traykovski  at  Atchafalaya  shelf  using  both  IDV  and  2D  models. 

IMPACT/APPLICATIONS 

The  present  research  efforts  focus  on  developing  a  detailed  numerical  modeling  framework  for 
cohesive  sediment  transport.  Currently,  the  Pis  are  also  actively  participating  in  other  related 
initiatives.  It  is  expected  that  by  the  end  of  this  2-year  project,  a  numerical  tool  will  be  made  available 
to  the  research  community  for  small-scale  cohesive  sediment  transport  processes.  The  Pis  are  also 
partners  of  the  ongoing  NOPP-CSTM.  In  this  collaborative  modeling  project,  we  utilize  the  small-scale 
model  to  study  wave-supported  gravity-driven  fluid  mud  and  provide  parameterizations  for  large-scale 
coastal  modeling  system. 
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Figure  1:  (a)  Snapshot  of  the  model  results  for  free-surface  elevation,  wave  velocity 
field  and  mud  concentration  in  the  numerical  wave  tank.  The  dashed  curve  near  the 
free-surface  is  the  corresponding  free-surface  elevation  without  fluid  mud.  (b)  Enlarge 
view  near  the  bed  with  mud  concentration,  wave  boundary  layer  velocity  at  three 
locations.  The  dashed  curves  are  corresponding  velocity  profiles  without  the  presence 
of  mud.  (c)  Comparison  of  wave  guage  data  atx=190  m.  Dashed  curve  represents 
model  results  without  mud.  The  light  dotted  curve  is  guage  results  at  x=10m  (before 
propagating  into  the  mud  regime).  At  this  instant  (t=72  sec),  the  wave  height  is  damped 
by  18%  due  to  propagating  over  fluid  mud  regime  of  180  m  in  length. 
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